Solar Radiation Pressure and Deviations from Keplerian Orbits 
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Newtonian gravity and general relativity give exactly the same expression for the period of an 
object in circular orbit around a static central mass. However, when the effects of the curvature 
of spacetime and solar radiation pressure are considered simultaneously for a solar sail propelled 
satellite, there is a deviation from Kepler's third law. It is shown that solar radiation pressure affects 
the period of this satellite in two ways: by effectively decreasing the solar mass, thereby increasing 
the period, and by enhancing the effects of other phenomena, potentially rendering some of them 
detectable. In particular, we consider deviations from Keplerian orbits due to spacetime curvature, 
frame dragging from the rotation of the sun, the oblateness of the sun, a possible net electric charge 
of the sun, and a very small positive cosmological constant. 
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In the last decade, the observation and analysis of 
satellite motion has provided an abundance of data with 
which to test basic physical principles. Examples include 
the Pioneer anomaly, which is an unexplained accelera- 
tion of the Pioneer 10 and 11 spacecraft on escape tra- 
jectories from the outer solar system 0, Q , and the flyby 
anomaly, for which the velocities of the Galileo, NEAR 
and Cassini spacecraft are different from what is expected 
after Earth flybys 0, H- In fact, the difficulties of ex- 
plaining these anomalies within the framework of stan- 
dard physics became a motivation to speculate on the 
unlikely possibility that they originate from new physics. 
Missions have been proposed @,1E 0| to further explore 
these anomalies, in order to better understand the laws 
of fundamental physics as they affect dynamics within 
the solar system. 

One of the most basic laws that describes motion in 
the solar system is Kepler's third law, which can be de- 
rived from Newton's law of gravitation and provides a 
relationship between the period T and the orbital radius 
r of an object orbiting the sun. Namely, T 2 oc r 3 with 
the proportionality constant given by 4tt 2 /GM, where G 
is the gravitational constant and M is the mass of the 
sun. Below, we discuss deviations from Keplerian orbits 
due to phenomena within conventional physics, some of 
which may be observed from the motion of solar sail pro- 
pelled (SSP) satellites 

In the Newtonian approximation, the sun is the source 
of a gravitational force on other masses. In the general 
relativistic framework, objects follow geodesies on the 
curved spacetime in the vicinity of the sun. At the same 
time, the sun is also a source of solar electromagnetic ra- 
diation, which produces an external force on objects via 
the solar radiation pressure (SRP). We will assume that 
the backreaction of the radiation on spacetime is negli- 
gible. Therefore, we can say that objects move in the 
photo-gravitational field of the Sun. 

The purpose of this Letter is to point out various 
sources of deviations from Keplerian orbits and discuss 



how the resulting change in period is enhanced by the 
SRP to the degree in which it may be observed for some 
cases. The phenomena discussed include the curvature 
of spacetime in the vicinity of the sun, described by the 
Schwarzschild metric, frame dragging due to the rota- 
tion of the sun, for which the curved spacetime is de- 
scribed approximately by the large-distance limit of the 
Kerr metric [9(, the oblateness of the sun, the effect of 
a possible small net electric charge on the sun, and a 
small positive cosmological constant. This list is by no 
means exhaustive. Some effects that we do not consider 
the gravitational field of planets, the magnetic field of 
the sun and the solar wind. While the reality is that 
all effects occur simultaneously, particular effects may be 
isolated by considering a variety of orbits. However, in 
order to get an approximate idea of the relative impor- 
tance of some effects in the presence of the SRP, we will 
consider them individually for the simple scenario of cir- 
cular orbits. 

According to Maxwell's electromagnetic theory, elec- 
tromagnetic waves carry the energy and linear momen- 
tum and the radiation pressure P exerted on a surface 
due to momentum transport by photons from the sun is 
given by 
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where S is the magnitude of the Poynting vector and the 
solar luminosity is L s = 3.842 x 10 26 W. Also, 0.5 < 17 < 
1, where r\ = 0.5 corresponds to the total absorption of 
photons by the satellite and 77 = 1 corresponds to total 
reflection. The resulting force is F = PA, where A is the 
surface area facing the sun. Thus, the acceleration due 
to this force can be expressed as 
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where m is the mass of the object. In fact, the mass per 
area a is a key design parameter for solar sails @, Gi| • 
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We will first consider the effect of the SRP on Keple- 
rian orbits in the Newtonian approximation for gravity. 
The SRP force is repulsive and the gravitational force 
is attractive. Both forces act along the same line (we 
are restricting ourselves to the case in which the sur- 
face is directly facing the sun) and they both fall off as 
1/r 2 . Therefore, the consideration of both forces leads 
to a modification of the effective mass of the sun in Ke- 
pler's third law. Namely, the mass of the sun, which is 
M = 1.99 x 10 30 kg, can be effectively renormalized as 
M = M — k/G. The modified Kepler's third law can be 
expressed as 



4tt 2 
GM 



(3) 



As a first example, we will consider Mercury, whose 
average orbital radius is r — 5.79 x 10 10 m, which cor- 
responds to a period of about 88 days. Mercury has a 
mass of 3.30 x 10 23 kg and a radius tm — 2.44 x 10 6 m. 
This yields a = 1.76 x 10 10 kg/m 2 , where the effective 
area is nr 2 M . If we assume that no sunlight is reflected 
by Mercury, so that r\ = 0.5, then the increase in period 
is on the order of 10 -7 s which, as to be expected for any 
planet, is negligible. 

We will now consider a conventional satellite orbiting 
the sun at a distance of r = 1 AU« 1.50 x 10 11 m, which 
corresponds to a period of one year. If the mass of the 
satellite is 1000 kg and its area is 2 m 2 , then a = 500 
kg/m 2 . If we suppose that rj = 0.75, then the increase in 
period due to SRP is about 36 s, which could be observed. 

In the remainder of the Letter, we will consider an SSP 
satellite with the following specifications: 

r = 0.05 AU w 7.48 x 10 9 m, 

i] = 0.85 , cr = 1.31 x 10~ 3 kg/m 2 . (4) 

If the acceleration due to the SRP is ignored, then the 
corresponding orbital period would be about 4 days. 
When the SRP is taken into account, the period is about 
70 days. 

We will now consider the simultaneous effects of the 
SRP and curved spacetime in the vicinity of the sun. 
We assume that the backreaction of the electromagnetic 
radiation on the background geometry is negligible so 
that it acts on the SSP satellite only via the SRP. The 
static exterior spacetime of the sun is described by the 
Schwarzschild metric, which is given by 



ds 2 = -fc 2 dt 2 + f- l dr 2 
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where 
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dVL 2 = d9 2 + sin 2 9 d(f> 2 . (6) 



Spherical symmetry allows us to orient the coordinate 
system so that the orbit is confined to the equatorial 
plane at 6 = tt/2, and thus pg — 0. Since the metric is 
independent of time and the azimuthal direction <fi, the 



corresponding components pt and p^ of the 4- momentum 
are conserved. We define the constants of motion E = 
—pt/m and L = p<j,/m, where m is the rest mass of the 
SSP satellite. Note that E and L are the energy and 
angular momentum of a unit mass, respectively. Thus, 
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where r is the proper time. In the absence of the SRP, 
p 2 = —m 2 c 2 yields 
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Differentiation of ([5]) with respect to r gives the radial 
component of the 4-acceleration 



a r-fl + - — + 3GML2 
We will now turn on the SRP, so that 
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where k is given in Q . Note that even though the coordi- 
nate r does not measure the proper distance, the surface 
area of a sphere is still given by 47rr 2 , which means that 
the magnitude of the Poynting vector as well as the ac- 
celeration are given by the same expressions as in the 
Newtonian approximation. Equating the expressions for 
a r given in ((9]) and (fTTJ|) and taking the first integral gives 
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Note that the SRP reduces the effective mass only in the 
term which is present for Newtonian gravity. 
For circular orbits, 



This yields 
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Using dt/dr — p t /m and dtfi/dr — p^/m, we find the 
orbital period T to be given by the expression 
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Keeping only the leading correction due to the curvature 
of spacetime gives 
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Notice that the simultaneous effects of the SRP and 
curved spacetime lead to a deviation from Kepler's third 
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law, which is completely absent without the SRP. For 
the specification given in we find that this yields an 
increase in the period of about 0.6 s. 

The rotation of the sun causes frame dragging, which 
affects the trajectories of orbiting objects. The external 
spacetime of a slowly rotating body with mass M and 
angular momentum J is described approximately by the 
large-distance limit of the Kerr metric [9J] 

AC J rlr 2 

ds 2 = -fc 2 dt 2 - sin 2 dtd(j> +^-+ r W, (16) 
c z r j 

where / is given by ((6|). We do not use the full Kerr 
metric since it does not seem to describe the external 
spacetime of a rotating material body, because it does 
not smoothly fit onto metrics which describe the interior 
region occupied by physical matter. Since there are cor- 
rections to this metric in higher-order J we will work up 
to only linear order in J, which suffices for the slowly 
rotating sun. Note that J > for a prograde orbit with 
respect to the sun, while J < for a retrograde orbit. 

Restricting ourselves to orbits that lie within the equa- 
torial plane, we define the constants of motion E = 
—p t /m and L = p^/m dUdl]. We have 



(17) 
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and p r = m dr/dr. In the absence of the SRP, p 2 
—c 2 m 2 yields 
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Differentiating this with respect to r gives 
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Turning on the SRP, the equation of motion is given by 
(TTUJ), which is not altered by linear terms in J. Equating 
the expressions for a r given in (fT0|) and (fT9|) and taking 
the first integral gives 
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v dr / r 
Applying the conditions for a circular orbit given in (|12p . 
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where X = c^r 2 - c 2 G{3M + M)r + AG 2 MM. Using 
dt/dr — p t /m and dcfr/dr — p^/m, we find the orbital 



period to be 
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0(J 2 ). 



Keep only the leading contributions due to spacetime cur- 
vature and frame dragging gives 
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The speed of the outer layer of the sun at its equator is 
v sa 2000 m/s at the equatorial radius R w 7 x 10 8 m. 
If we assume that the core of the sun rotates with the 
same angular speed, then J = ^MvR ps 10 42 kg m 2 /s. 
Without the effects of the SRP, frame dragging leads to 
an increase (decrease) in the period of only 4 x 10~ 5 s 
for a prograde (retrograde) orbit. With the SRP and the 
specifications in (U) , frame dragging leads to a change in 
the period of about 0.01 s. 

We will now consider the effect of the oblateness of the 
sun. Working in Newtonian gravity, the external gravi- 
tational potential of an oblate spheroid is given by [l3T | 
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where J n are the multipole mass moments, R is the equa- 
torial radius of the sun and P n are the Legendre polyno- 
mials. Note that the effective mass in the first term is 
renormalized by the solar radiation pressure, whereas the 
multipole mass moments are not affected. We will con- 
sider only the case for which the orbit of the SSP satellite 
is confined to the equatorial plane. Then the acceleration 
is purely in the radial direction, given by 
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This leads to the relation 
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where we consider the first sub-leading term due to the 
oblateness. The oblateness of the sun has recently been 
measured with unprecedented precision to be as much as 
J 2 w 9x 10~ 6 during active phases of the solar cycle [l4T |. 
Without the SRP, the oblateness of the sun increases the 
period by about 0.02 s. With the SRP and the speci- 
fications given in the sun's oblateness results in an 
increase in the period of about 105 s. It can be shown 
that the next mass moment J4 increases the period by 
only about 0.0006 s. 

We will now consider the effect of a small amount of 
net charge Q. It has been suggested that the sun has 
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a net charge of up to Q w 77 C [15|. The spacetime is 
described by the Reissner-Nordstrom metric, which has 
the form {5]), where the function / is now given by 



/ = 1 
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and fc is the Coulomb constant. Taking into account the 
electric force, we find that 
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where q is the charge of the solar sail. Note the Q 2 term 
in the period, due to the backreaction from the sun's 
charge on the geometry, is present even for a neutral SSP 
satellite. However, this term increases the period by an 
amount of only 1CT 35 s (10 -38 s without the SRP), which 
reflects the fact that the backreaction of the charge on 
the geometry is negligible, as to be expected. 

Primarily due to the photoelectric effect, but also the 
Compton effect and electron-positron pair production, a 
solar sail made from Beryllium, for example, will equili- 
brate to a charge per area of 0.065 C/m 2 [H, If the 
SSP satellite has a mass of 1000 kg, this gives a charge 
q = 5 x 10 4 C and an increase in period of about 230 s 
(0.05 s without the SRP). Thus, due to the SRP, even a 
small net charge Q could yield a measurable increase in 
the period. However, the challenge would be to isolate 
the effects of the sun's highly variable magnetic field. 

Supernovae observations suggest that our universe 
mig ht have a very small positive cosmological constant 
[l8l [l9j . Since the SRP enhances a variety of small ef- 
fects, one could ask how much the effects of a cosmolog- 
ical constant are enhanced. In the presence of a cosmo- 
logical constant A, the spacetime in the vicinity of the 
sun is described by a metric of the form ([5]), where the 



function / is now given by 
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Considering the effect of A for r ^> GM/c 2 , we get 
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For A « 10 m , this leads to an insignificant increase 
in period of 10~ 17 s (10~ 21 s without the SRP). While 
one might be able specify values of a and r such that 
the correction factor is actually significant, the period 
would then be too large to make this a feasible test for 
the presence of a cosmological constant. 

Below we summarize our results for the change in pe- 
riod AT of an SSP satellite with specifications j4} . 



Phenomenon 


AT without SRP 


AT with SRP 


Spacetime Curvature 


s 


0.6 s 


Frame Dragging 


±4 x 10~ 5 s 


±0.01 s 


Oblateness of Sun 


0.02 s 


105 s 


Net Charge of Sun 


0.05 s 


230 s 


Cosmological Const. 


10- 21 s 


10- 17 s 



The effects of some phenomena are increased by sev- 
eral orders of magnitude due to the SRP. While this 
could be sufficient to observe certain phenomena, the 
SRP would also enhance additional effects that we have 
not discussed, and a variety of orbits would be needed 
in order to isolate particular effects. Observations of the 
above deviations in the period would provide an interest- 
ing confirmation of these phenomena. 
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